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Abstract 

Once super symmetric neutralinos x° are produced copiously at e + e~ linear collid- 
r> \ ers, their characteristics can be measured with high precision. In particular, the 

& ' fundamental parameters in the gaugino/higgsino sector of the minimal supersym- 

metric extension of the standard model (MSSM) can be analyzed. Here we focus on 
the determination of possible CP-odd phases of these parameters. To that end, we 
exploit the electron/positron beam polarization, including transverse polarization, 
as well as the spin/angular correlations of the neutralino production e + e~ — > XiXj 
and subsequent 2-body decays x° - > X^h, x\^i &r^i using (partly) optimized CP- 
odd observables. If no final-state polarizations are measured, the Z and h modes 
are independent of the Xi polarization, but CP-odd observables constructed from 
the leptonic decay mode can help in reconstructing the neutralino sector of the 
CP-noninvariant MSSM. In this situation, transverse beam polarization does not 
seem to be particularly useful in probing explicit CP violation in the neutralino 
sector of the MSSM. This can most easily be accomplished using longitudinal beam 
polarization. 
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1 Introduction 



In the minimal supersymmetric standard model (MSSM) [1], the spin-1/2 partners of the 
neutral gauge bosons, B and W3, and of the neutral Higgs bosons, H® and H®, mix to 
form the neutralino mass eigenstates Xi (^ = 1;2,3,4). The corresponding mass matrix in 
the (B,W 3 ,H?, Hi) basis 
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(1) 



contains several fundamental supersymmetry parameters: the U(l) and SU(2) gaugino 
masses Mi and M 2 , the higgsino mass parameter /x, and the ratio tan/5 = V2/V1 of the 
vacuum expectation values of the two neutral Higgs fields. Here, sp = sin/3, cp = cos/3 
and Sw,cw are the sine and cosine of the electroweak mixing angle 6w 

In CP-noninvariant theories, the mass parameters Mi )2 and /x are complex. By re- 
parameterizing the fields, M 2 can be taken real and positive without loss of generality. 
Two remaining non-trivial phases are attributed to Mi and /x: 



Mi = I Mi I e 



and 



/x — |/x| e 



(0<$i,$ m <2tt). 



(2) 



The existence of CP-violating phases in supersymmetric theories induces, in general, 
electric dipole moments (EDM) [2]. The current experimental bounds on the EDM's 
constrain the parameter space including many parameters outside the neutralino/ chargino 
sector [3]. Detailed analyses of the electron EDM show [3, 4] that the phase $ M must 
be quite small, unless selectrons are very heavy* In contrast, large values of $1 are 
allowed even for rather small selectron masses. The CP-violating phase $1 can therefore 
play a significant role in the production and decay of neutralinos, which is most easily 
investigated at (linear) e + e~ colliders [5, 6, 7, 4, 8]. 

Neutralinos are produced in e + e~ collisions, either in diagonal or mixed pairs [9]. If 
the collider energy is high enough to produce all four neutralino states, the underlying 
SUSY parameters {|Mi|, $1, M 2 , \fj,\, tan/3} can be extracted from the masses m^o 
(z=l, 2,3,4) and the cross sections [10, 11]. At the first stage of operations of a linear e + e~ 
collider, however, only the lighter neutralinos may be accessible. If X1X2 * s the only visible 
neutralino pair that is accessible, measuring their masses and (polarized) production 
cross sections may not suffice to determine the parameters of the neutralino mass matrix 
completely; the detailed analysis of x° decays will then be very useful. Moreover, even 
if sufficiently many different states are accessible to determine all the parameters 
appearing in Eq. (1), analyses of neutralino decay will offer valuable redundancy. After 

Targe values of can also be tolerated for moderate selectron masses if tan (3 is close to 1. However, 
this possibility is essentially excluded by Higgs boson searches at LEP. 
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Figure 1: Feynman diagrams for five mechanisms contributing to the production of di- 
agonal and non-diagonal neutralino pairs in e + e~ annihilation, e + e~ — > XiXj (hj = l~4)- 

all, a theory can only be said to be tested successfully if experiments over-constrain its 
parameters. 

In the present work we systematically investigate, both analytically and numerically, 
the usefulness of electron and positron beam polarization, including transverse polariza- 
tion, for the analysis of neutralino production and decay at e + e~ colliders. To this end, 
we exploit spin/angular correlations of the neutralino production e + e~ — > X2X1 an d sub- 
sequent two-body decays of X2 — > X±h, and X2 — * ^ ± ^ T followed by £ ± — > i^Xi for 
probing the CP properties of the neutralino sector in the MSSM. Due to the Majorana 
nature of neutralinos, the decay distributions of two-body decays X2 ~^ Xi^, Xi% are 
independent of the X2 polarization, unless the polarization of the Z boson is measured. 
These modes can still be used to probe a production-level CP-odd asymmetry, which 
however turns out to be small in the MSSM. The slepton mode X2 ~^ * s an opti- 
mal polarization analyzer of the decaying neutralino. We can construct several CP-odd 
"decay" asymmetries that are sensitive to the X2 polarization vector. Our main emphasis 
is on observables that fully reflect the non-trivial angular dependence of CP-odd terms, 
except for the angular dependence appearing in the propagators. Although they are not 
perfectly optimal, these CP-odd asymmetries have much higher statistical significance 
than the conventional ones, as demonstrated with numerical examples below. 

The remainder of this article is organized as follows. Section 2 describes neutralino 
production, including the polarization of the neutralinos, for arbitrary beam polarization. 
Two-body decays of polarized neutralinos are discussed in Sec. 3. Section 4 deals with 
the reconstruction of X1X2 final states with invisible Xi- The formalism of "effective 
asymmetries" is described in Sec. 5, and numerical examples for these asymmetries are 
shown in Sec. 6. Finally, Section 7 contains a brief summary and some conclusions. 

2 Neutralino production in e + e~ collisions 

The neutralino pair production processes in e + e~ collisions 

e-(p, a) + e + (p, a) -> # (Pi, \) + X°(Pj, ^) ( i, j = 1, 2, 3, 4) (3) 

2 



are generated by the five mechanisms of the Feynman diagrams in Fig. 1, with s-channel 
Z exchange, or t- or w-channel &l,r exchange. Here a, a, Aj, and Xj denote helicities. For 
the analytical calculation, we take a coordinate system where the production occurs in the 
(x, z) plane and the incident electron beam moves into +z direction. The four-momenta 
appearing in Eq. (3) are then given by 

p = ^(1,0,0, 1), 



Pi 
Pj 



^(1, 
2 v ' 

^-(e h A 1/2 sin6, 0, A 1/2 cos6) 



p = ^-(1,0,0,-1) 
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^(e,-, -A 1 / 2 sine, 0, -A 1/2 cos 6) , (4) 



where 



e, = 1 + $ - $ , ej = 1 + /i 2 - , 

to* = m x% 3 l^~ s > A = i 1 ~ f*i ~ v)) 2 ~ ^ifij ■ ( 5 ) 



The transition matrix element, after an appropriate Fierz transformation of the cl,r 
exchange amplitudes, can be expressed in terms of four generalized bilinear charges Q a p- 



e 2 



T (e + e~ - = ~ Q*P [v(e + h,P a u(e~)] [uffiWPfivtfj)] ■ (6) 

These generalized charges correspond to independent helicity amplitudes which describe 
the neutralino production processes for completely (longitudinally) polarized electrons 
and positrons, neglecting the electron mass as well as eL-e# mixing.* They are defined in 
terms of the lepton and neutralino couplings as well as the propagators of the exchanged 
(s)particles [6, 11]: 

s w c w 
D z 

QrL = +^2~ + DtR9Rij, 



c w 
D 



QLR 2—2- (S\V - \ + DtLglij, 

s w c w 
Dz 

Qrr = — 2~Zij - D U Rg Rij . (7) 
c w 



* fL~fn mixing is proportional torn/ unless one tolerates deeper minima of the scalar potential where 
charged sfermion fields obtain nonvanishing vacuum expectation values; although it can be enhanced at 
large tan (3 or for large trilinear A— parameters, selectron mixing is generally negligible for collider physics 
purposes. 
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The first index in Q a p refers to the chirality of the current, the second index to the 
chirality of the x° current. The first term in each bilinear charge is generated by Z- 
exchange and the second term by selectron exchange; Dz, D t L,R and D u l,r respectively 
denote the s-channel Z propagator and the t- and it-channel left /right-type selectron 
propagators: 

Z s — m\ + irrizTz ' 
g 

D t L,R = ~ o aIld t^U, (8) 

t — mi 

with s — (p + p) 2 , t = (p — pi) 2 and u = {p — pj) 2 . The matrices Zij, guj and gmj can be 
computed from the matrix N diagonalizing the neutralino mass matrix [1] 

^• = (awv; 3 -awv; 4 )/2, 

guj = (N i2 c w + NnSw)(N* 2 c w + N^sw) / As^Cw , 
gRij^NaNJJcly. (9) 
They satisfy the hermiticity relations of 

Zij = Zji > 9uj = glji > gmj - g*Rji ■ 

If the decay width V z is neglected in the Z boson propagator Dz, the bilinear charges 
Q a/ 3 satisfy similar relations, Qapix^, X°j, t, u ) — Qapix^ Xi > u i t)- These relations are very 
useful in classifying CP-even and CP-odd observables. 



2.1 Production helicity amplitudes 

With the e 1 * 1 mass neglected, the matrix element in Eq. (6) is nonzero only if the electron 
helicity is opposite to the positron helicity. We write the helicity amplitudes as 

T(<7, a, Aj, \j) = T{a, -a, A i7 \j) S & _ CT = 2na (a; A» Xj) S & _ CT , (11) 

where a, Aj, Aj = ±. Explicit expressions for these helicity amplitudes are [6]: 

(+; ++) = - [Q RRy /Vi+Vj- + QRLVVi-Vj+] sin , 

(+; H — ) = - [Q R Ry/r]i+'nj+ + QRLy/Vi-Vj-] (l + cos @) , 

( + ! -+) = + [QRRy/Vi-Vj~ + QRL^Vi+Vj+] (1 - COS @) > 

(+; ) = + [QRRy/Vi-Vj+ + QRL^/Vi+Vj-] sin 6 , 

(-; ++) = - [QLLy/Vi-Vj+ + QiR^Vi+Vj-] sin 6 , 

(-; +-) = + [QLLy/Vi-Vj- + QLRy/V l +Vj+] i 1 ~ COS 6) , 
(-; - + ) = - [QLLy/Vi+Vj+ + QLR^/Vi-Vj-] (1 + COS 6) , 

(-; ) = + [QLLy/Vi+Vj- + QlRy/Vi-Vj+l SU1 6 , (12) 

where r]i± = e$ ± A 1 / 2 and rjj± = ej ± A 1 / 2 . In the high energy asymptotic limit, r]i + and 
r\i- approach 1 and 0, respectively; only the helicity amplitudes with opposite x° and x° 
helicities survive. 
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2.2 Production cross sections 

We analyze neutralino production for general e^ 1 polarization states. With the scattering 
plane fixed as the (x, z) plane, the azimuthal scattering angle appears in the description 
of the e polarization vectors: 

P* e - = (P T cos$,-P T sin$,P L ), ~P e+ = (P r cos(77-$),P T sin(77-$),-P L ), (13) 

where rj is the relative angle between the transverse components of two polarization vec- 
tors. The density matrices p (p) of the electron (positron) in the {+, — } helicity basis are 
[13] 

_1/1 + P L P T e^ \ __lf 1 + P L -P T e-^-*) \ 
P "2 VPre-** 1-Pl)' P ~ 2 { -P T e^ 1 -P L )' [ } 

The polarized differential cross section is given by 

da A 1 / 



dQ 64vr 2 s 



T\\ (15) 



where 



|T| 2 = T(a,(T,X i ,Xj)T*(a / ,a , ,Xi,Xj)p mTl p 9 , 9 . (16) 

Note that the order of indices of TVs is opposite of that of p CTO -/ due to the difference 
between the particle and the antiparticle. Inserting Eqs. (12) and (14) into Eq. (16) 
yields 

%m = Al/2 - PlPl ^ ^ + - ^ 



+ P r P T cos(2$-77)Eg T + P r P r sin(2$-77)Eg JV , (17) 



where 



Egc = [1 - {pi - p 2 ) 2 + A cos 2 6] Q x + 4 W Q 2 + 2A 1 / 2 Q 3 cos 0, 

sg £ = [i - (/i 2 - /i 2 ) 2 + a cos 2 e] q; + 4 Atj/ i j g / 2 + 2\ 1 ' 2 q' z cos e, 

T>u T = A Q 5 sin 2 9, 

S^ = -Ag^sin 2 6. (18) 

Expressions for all relevant quartic charges in terms of bilinear charges Q a p are given 
in Table 1, which also lists the transformation properties under P and CP. Non-zero 
transverse beam polarization allows to probe four new quartic charges, Q$, Qq, Q' 5 , 
and Q' 6 . 
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Table 1: The independent quartic charges describing e + e — > XiXj- 



T 1 

P 


CP 


Quartic charges 


even 


even 


D — 
VI — 




\n r j 2 i 






D — 

V2 — 


gJte [^RR<^ RL + ^LL^LRl 








Qs = 


\\\Qrr\ 2 + \Qll? -\Qrl\ 2 ~ 

A L 1 Ail | 1 | ^/ -L/.L/ \ \ ^< 1X.1L/ \ 


IQii?| 2 l 

| J-/-IX | J 






g 5 = 


\Ue [QrrQ* lr + QllQr L ] 






odd. 


Qa = 


o^ym fQRRQRr + QllQ*t r1 








Qe = 






odd 


even 


Q'i = 


1[|Qh«I 2 + IQhl| 2 -|Qll| 2 - 


IQl^I 2 ] 






Q'2 = 


2^ e [QrrQrl - QllQIr] 








Qs = 


\[\Qrr\ 2 + \Qlr\ 2 -\Qll\ 2 - 


\Qrl\ 2 ] 






Q' 5 = 


2^e [QrrQ* L r - QllQrl] 






odd 


Q',= 


l^m [QrrQ* rl — QllQlr] 








<% = 


i^m [QrrQ* lr — QllQrl] 





2.3 Neutralino polarization vector 

The polarization vector V 1 = {V?, V l N , V L ) of the neutralino Xi is defined in its rest frame. 
The longitudinal component V\ is parallel to the xl flight direction in the cm. frame, V % T 
is in the production plane, and V l N is normal to the production plane. In order to extract 
the vector V 1 , we first define the polarization density matrix for the out-going neutralino 
y°- 

Explicit expressions for the helicity amplitudes (a; A,/Aj) are given in Eq. (12). The po- 
larization vector of the neutralino x® is then given by 

P = Tr( oV) = -L (A?, A« , A?) . (20) 
A u 

We can decompose the three polarization components as well as the unpolarized part 
according to combinations of e ± polarizations: 

A^ = (1 - P L Pl)^uu + ( P L - Pl)^ul + PtPt{^ut C (2^-v) + ^UN S (2^>-7 1 )} , 
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Ag = (1 - PlPlMu + (Pl ~ Pl)^L + PtPt{AtC{2*-v) + , 

A? = (1 - P L P L )^4u +(Pl~ Pl)^l + PriM^c^^) + E^ N s {2 ^ v) } , 

Ag = (l-P L P L )Z% u + (P L -P L )^ L + P T P T {Z% T c^ rl) + (21) 

where c^.^ = cos(2$ — rj), S( 2 $_^) = sin(2$ — rj), and the (£> = U, L, T, N) are in 
Eq. (18). The Y^bu, which survive even without beam polarization, are given by 

T^Ijj = 2(1 - /x 2 - nfj cos 9 Qi + Apuij cos G Q 2 + A 1/2 {1 + cos 2 6 - sin 2 6( / u 2 - /i 2 )} Q' 3 , 

= -2 sin 6 [{(1 - ^ 2 + /i 2 ) Q[ + A 1 / 2 cos 6 Q'^pa + (1 + ^ 2 - Q' 2 ] , 

= 2A 1/2 /x j sin6Q 4 - (22) 

The remaining £ab, which contribute only with non-trivial e* 1 polarization, are 

^ll — [A + 1 — (yw? — /i 2 ) 2 ] cos 6 Qi + 4/ii/ij cos 9 Q 2 
+A 1/2 [1 + cos 2 6 - sin 2 6 (/i 2 - n))} Q 3 , 
S^ T = A^(l + ^-^) S in 2 0Q' 5 , 



^ = -A 1/2 (l + /i 2 -^ 2 )sin 2 eg 
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E? £ = -2 sin 6 { [(1 - /i 2 + /i 2 ) Qi + A 1/2 cos 9 Q 3 K + (1 + /x? - Qi} , 
Y% T = A 1/2 /iiSin2eg' 5 , 

= -A 1/2 /i iS in20g 6 , 
S^ L = 2A 1 / 2 /!,- sinGQl, 
^nt - -2A 1/2 ;Ui sin0Q 6 , 

= -2A 1 / 2 //* sin0Q' 5 . (23) 

The P and CP properties of all these quantities are identical to those of the quartic charges 
in Table 1. In particular, the five quantities T, LN , T, TN , T>nu an d Svl are CP-odd. 

Brief comments on the reference frame are in order here. In the coordinate system 
which we have employed so far, the scattering plane is fixed, while the direction of e ± 
transverse polarization vectors differs from event to event. For a real experiment, fixed 

polarization vectors should be more convenient. We define the transverse part of P e - 
as +x direction; the x and y components of the outgoing neutralino four-momentum 
Pi are then proportional to cos $ and sin $, respectively. In this coordinate system the 
scattering plane changes from event to event. Since only the relative angles between the 

polarization vectors and the scattering plane are relevant, the final results in Eqs. (17) 
and (21) are still valid. In this new coordinate frame, the xl polarization vector can be 
explicitly written as 

V 1 = V T e T + V N e N + V l L e L , (24) 

where the following three unit vectors form a co-moving orthonormal basis of the three- 
dimensional space: 

tT = (cos $ cos 0, sin $ cos 9, — sin 0) , 
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e/v = (— sin $, cos $, 0) , 

ex, = (cos $ sin 6, sin $ sin 6, cos 0) . (25) 

Probing CP violation in the MSSM neutralino sector involves the four quartic charges 
Qi, Q'4, Qq and Q' 6 for i 7^ j. Their characteristic features can be analytically understood 
from their explicit expressions in terms of the neutralino mixing matrix N. With Yz 
neglected in the high energy limit, they are 



<ff = ^iV K =F (*w ~ V2) 2 ] 
D 



>j ' 



2r 2 



s 2 - 1/2 

D tR + D uR )Sm(Z i:j g Rij ) ± ~ 2 (D tL + D UL )^m(Z ij g Lij ) 



s w 



+-D uR D tR $m(g 2 Rij ) T ^D uL D tL %m{g 2 Lij ) , 

Q { 6 = ^rD z (D tL ± D uL )Zm(Z ljg * Uj ) + S \ ~ \ /2 D z (D uR ± D tR )^m(Z ij9 * Rij ) 

+ l -(D uR D tL ± D tR D uL )$m( gLij g* Rij ) , (26) 

where the explicit form of Zij, guj and g R jj are listed in Eq. (9). From the propagator 
combinations, we see that the quartic charge Q' & is forward-backward asymmetric with 
respect to the scattering angle while the other three quartic charges, Q4 and Qq, are 
forward-backward symmetric. 

The relative sizes of the four CP-violating quartic charges indicate which observables 
should be promising to investigate experimentally. Let us first consider the generic case of 
small gaugino-higgsino mixing (with substantial CP phase $1). Small mixing is generally 
obtained if the entries in the off-diagonal 2x2 blocks in the neutralino mass matrix 
are smaller than those in the diagonal blocks, allowing an expansion in powers of mz- 
Analytic expressions for N using this expansion, given in Ref. [4], help to estimate the 
sizes of the Q Q%\. In particular, the last term contributing to Q4 in Eq. (26), which is 
proportional to sin$i, is not suppressed by small mixing angles: Q4 and Q' 4 survive even 
without any gaugino-higgsino mixing. In contrast Qq and Q' 6 only start at 0{m 2 z ). This is 
related to the observation that, in the notation of Ref. [11], Q$ probe Dirac-type phases, 
which vanish in the absence of nontrivial mixing between neutralino current eigenstates, 
whereas Q 4 probe Majorana-type phases, which survive in this limit. In the generic case 
of small gaugino-higgsino mixing, therefore, the size of Q 4 is much larger than that of 
Qf. In the case of strong gaugino-higgsino mixing, however, Q % , which can only be 
probed with transversely polarized beams, could exceed Q4 and/or Q' A . 
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3 Two— body neutralino decays 



The decay patterns of heavy neutralinos (Xi>x) depend on their masses and the masses 
and couplings of other sparticles and Higgs bosons. In this article we focus on the two- 
body decays of neutralinos. It is possible that the kinematics prohibits some two-body 
tree-level decays. However, a sufficiently heavy neutralino can decay via tree-level two- 
body channels containing a Z or a Higgs boson and a lighter neutralino [14], and/or into 
a sfermion-matter fermion pair. 

Of particular interest in the present work are the following two-body decay modes: 
X°^X° k Z, Xi^xth and Xi - , (27) 

with i = e or /x. If any of these processes is kinematically allowed, it will dominate any 
tree-level three-body decay. 

The relevant couplings are 

(r L \t R \x°) = +<^I4JX°>* = -V2<7%iV*, (£i\£i\x-) =0, (28) 
(x° kR \Z\x iR ) = -{xIl\Z\x'lY = [N i3 N* k3 - NhNM , 



2c w 



(x° kL \h&) = +(xl R \h\Xi L )* = § [(#w " t w N kl )(s a N i3 + c a N z ,) + (i <- k)] , 

where s a = cos a, c a = sin a, and a being the mixing angle between the two CP-even 
Higgs states in the MSSM [1]. Note that the Z coupling is proportional to the higgsino 
components of both participating neutralinos, whereas the Higgs coupling requires a hig- 
gsino component of one neutralino and a gaugino component of the other.* Since the 
lighter neutralino states Xi 2 are often gaugino-like, this pattern of couplings implies that 
Xi ~^ X\h decays will often dominate over the (kinematically preferred) xl ~~ * X\Z decays. 
However, the xi ~^ decays only depend on the gaugino components of the decaying 
neutralino. If kinematically accessible, they can have the largest branching ratios. 

Note also that the Majorana nature of neutralinos relates the left- and right-handed 
couplings of the Z and h boson to a neutralino pair; they are complex conjugate to each 
other, having an identical absolute magnitude. These relations lead to a characteristic 
property of the corresponding two-body decays, xl X k Z and Xi - > X k h: the decay 
distributions are independent of the polarization of the decaying neutralino x®, unless the 
polarization of the Z boson or x k is measured. In contrast, the slepton mode in Eq. (27) 
can be exploited as optimal polarization analyzer of the decaying neutralino, if the small 
lepton mass is ignored; as noted earlier, this implies that £l~£r mixing is ignored as well.^ 



*If Sm^ = ra^o — m^o 3> rnz, the decay into longitudinally polarized Z bosons gets enhanced by a 
factor (5m^/mz) 2 ■ If Sm^ ~ 0{mz), three-body decays X2 ~~ * Xiff m ay dominate over X2 X\Z 
decays if ^ m^; this does not happen in models where the entire sparticle spectrum is described by 
a small number of parameters. 

^Xi — > r-j^T^ decays, where fx— tr mixing can be important, have been analyzed in Refs. [7]. 
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Furthermore, the decay distributions are completely determined by the relevant par- 
ticle masses, as well as by the Xi polarization vector (in case of Xi -> £r£ t decay). More 
explicitly, the angular distribution in the rest frame of the decaying neutralino xl is 

(29) 



T x dQ* 4vr 

where £z,h = for the Z and h decay modes, and &± = =)=1 for Xi ~^ ^r^ t with k\ being 
the unit vector in p direction. The former two decay modes can probe only "production" 
asymmetries, whereas the (s)leptonic decay mode can probe "decay" asymmetries also, 
which are sensitive to the Xi polarization. 



4 Event reconstruction 

We focus on e + e~ — > X2X1 production, and assume Xi to be stable (or possibly to decay 
invisibly). The only visible final state particles therefore result from x 2 decay, which 
simplifies the analysis. Moreover, this is the kinematically most accessible neutralino pair 
production with visible final state; indeed, it is often the first sparticle production channel 
accessible at e + e~ colliders [15]. 

An important difference between x° ~^ Xi(h, Z) and x% ~ > ^r^ t ~ * Xi^ + ^~ is the de- 
gree of event reconstruction. The latter decay chain allows complete event reconstruction 
(with an, at least, two-fold ambiguity), whereas the former does not. This can be seen by 
counting unknowns. The XiXi(hj final states contain six unknown components of Xi 
momenta (we are assuming that the masses of all produced particles have already been 
determined [10], so that the energies can be computed from three-momenta); this has 
to be compared with four constraints from energy-momentum conservation, and a single 
mass constraint, (p^o + P(h,z)) 2 — m |o- One quantity remains undetermined. 

In contrast, XiXi^ + ^~ final states produced from an on-shell £^ have two invariant 
mass constraints. With an equal number of constraints and unknowns, the event can 
be reconstructed [8]. An explicit reconstruction may proceed as follows. Let k\ and k 2 
be the four-momenta of the two charged leptons in the final state, and pi and q the 
four-momenta of the two neutralinos; here k 2 and q originate from Ir decay. Note that 
the energy p\ is fixed from two-body kinematics, see Eq. (4). Then q° is determined 
from energy conservation, once the lepton energies are measured. The invariant mass 
constraint (k 2 + q) 2 = m~ can fix the scalar product k 2 ■ q. The second mass constraint 

(k\ + k 2 + q) 2 = m-o is used for k\ ■ q. When writing the unknown three-momentum q 

as q = ak\ + hk 2 + c(k\ x k 2 ), the two coefficients a and b can be computed from the two 
scalar products k 2 ■ q and k\ ■ q determined above; note that the term proportional to c 
drops out here. The last coefficient c can be computed from the known energy q° with 
two-fold ambiguity. 
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Once q is known, pi follows immediately from momentum conservation. We can read 
off the production angles G and $. This also allows to compute the X2 three-momentum 
f>2 = ki + k 2 + q = —pi (in the cm. frame). With the known x° energy, we boost into 
the X2 res t frame, and read off the X2 decay angles Q* and $*; recall that there is a 
non-trivial dependence on these decay angles via Eq. (29). 

So far we have assumed that we know which of the two charged leptons in the final 
state originates from the X2 decay, and which one from £r decay. Since, owing to its 
Majorana nature, x® wm decay into both £r£~ and £r£ + final states with equal branching 
ratios, the charge of the leptons does not help this discrimination of the origin of two 
charged leptons. A unique assignment is nevertheless possible if the two mass differences 
52R = m^o —ra- lR and Sri = mi R —tti^o are very different from each other: if 5 2 _r ^> Sri, the 
more energetic (harder) lepton will originate from the first step of X2 decay, and the less 
energetic (softer) lepton comes from £ R decay; if 5 2 r <C Sri the opposite assignment holds. 
However, if S2R — Sri, both assignments often lead to physical solutions if the procedure 
for event reconstruction outlined above is applied. In this unfavorable situation there is 
a four-fold ambiguity in the event reconstruction. 

Finally, we note that background events can be also reconstructed, in some cases again 
with two-fold ambiguity. The main backgrounds to X2 ~> Xi(Z,h) decays are e + e~ — > 
ZZ, Zh production with one Z decaying invisibly. The e + e~ — > ZZ{-^ vv£ + £~), W + W~(- 
£ + ui£~U{), £ + £~(— > £ + £~XiXi) are the main backgrounds to X1X2 ~^ £ + £~XiXi produc- 
tion.* We can obtain a pure sample of signal events by discarding all events that can be 
reconstructed as one of the background processes. This ignores the effects of measure- 
ment errors, beam energy spread (partly due to bremsstrahlung) , as well as initial state 
radiation, but should nevertheless give a reasonable indication of the effects of cuts that 
have to be imposed to isolate the signal. 



5 Effective asymmetries 

We are interested in constructing CP-odd observables. Schematically, they are written as 

f = j dn^f(n)xc, (30) 

where da/dQ is the differential cross section, C = j Ldt is the total integrated lumi- 
nosity, and f{Q) is a dimensionless function of phase space observables. Introducing the 
luminosity in Eq. (30) simplifies the statistical analysis as presented below. 

Simple asymmetries are constructed from the choice / = ±1, where the phase space 
region giving / = +1 is the CP-conjugate of that giving / = — 1 [5, 8]. While very 
straightforward, this choice usually does not yield the highest statistical significance. We 

*Note that we include supersymmetric slepton production as background, since it does not contribute 
to the CP-odd asymmetries we wish to analyze here. 



11 



decompose the differential cross section into CP-even and CP-odd terms: 



^ = E e ^ (e) ( Q )+E o ^ 0) ( fi )> ( 31 ) 

» 3 

where the and Oj are constant coefficients (products of couplings and possibly masses) 
while the and p ' are CP-even and CP-odd functions, respectively, of phase space 
variables. The optimal variable to extract the coefficient Oj is then proportional to 
[16]. 

In our case this would lead to very complicated observables, due to the non-trivial 
angular dependence of the selectron propagators Dn u v LjR \ in Eq. (7). Moreover, the 
optimal variables would depend on both selectron masses. For simplicity, we construct 
our CP-odd observables by fully including the angular dependence in the numerators of 
Eqs. (17), (18), (21), (22), (23) and (29), but ignoring the angular dependence in the 
propagators. 

For dimensionless /, the quantity F in Eq. (30) is also dimensionless. The statistical 
uncertainty of F is then given by 

a\F) = Cx J dfl^/ a (0). (32) 

This can be seen from the fact that C(da/dQ)dQ is the number of events in the phase 
space interval dQ. For the simple case of / = ±1, cr 2 (F) is simply the total number of 
events. With the quantity F and its statistical uncertainty cr(F), we can construct an 
effective asymmetry: 

A ^^Tc (33) 

Note that A is by construction independent of the luminosity. It is also invariant under 
transformations /(fi) — > cf(Q) for constant c, making A independent of the normalization 
of /. The statistical significance for A[f] is simply given by A[f] ■ 



6 Numerical analysis 

We are now ready to present some numerical results. We will first briefly discuss the 
relevant quartic charges that encode CP violation, before discussing "production" and 
"decay" asymmetries. 



6.1 Quartic charges 

Table 1 shows that the four quartic charges Q^, Qq, Q f 4 and Q' e are CP-odd. Equation (18) 
shows that Q' 6 is responsible for the production-level asymmetry, which requires transverse 
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beam polarization.* The remaining three CP-odd quartic charges can be probed only 
via the \\ polarization. Equations (22) and (23) show that Q4 contributes even for 
unpolarized e ± beams, whereas Q' A (Qe) only contributes in the presence of longitudinal 
(transverse) beam polarization. 
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Figure 2: The ratios of quartic charges Q4/Q1 (dotted green), Q'jQi (dashed blue), 
Qe/Qi (solid red) and Q'JQi (dot-dashed black). We fixed \M X \ = 0.5M 2 = 150 GeV, 
tan/3 = 5, mg L = 500 GeV and $ M = 0; the values of the other relevant parameters are 
as indicated in the figures. 

Figure 2 presents these four charges normalized to Qi, which largely determines the 
size of the unpolarized cross section far above threshold. All these ratios lie between — 1 

*We note in passing that the corresponding asymmetry for chargino production vanishes [17]: there 
is no equivalent of the e# exchange diagram, and the relevant 2x2 matrix diagonalizing the chargino 
mass matrix does not contain a reparametrization invariant phase. 
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and 1. We took \M X \ = 150 GeV, M 2 = 300 GeV (so that |Mi| and M 2 unify at the 
scale of Grand Unification [1]), a moderate tan/3 = 5, mg L = 500 GeV, and $ M = (as 
indicated by constraints on the electric dipole moments of the electron and neutron [2, 3]). 
The default choices of the other relevant parameters are \fi\ = 325 GeV, m e - R = 300 GeV, 
$i = 0.67T and y/s = 2E^ eam = 500 GeV, but one of these parameters is varied in each of 
the four frames of Fig. 2. Finally, we chose scattering angle cosG = l/v2 ; note that Q' 6 
vanishes at cos O = 0. 

The behavior of the curves in Fig. 2 can be understood with the help of the expressions 
in Eq. (26). The top-left frame shows the dependence of the four ratios on the phase $i. 
We see the typical behavior of CP-odd quantities, changing sign when sin$! changes 
sign, although not simple sine functions. Since we took to be close to M 2 , x° I s 
a strongly mixed state. However, xl is still mostly gaugino-like, so that |Z 12 | is quite 
small. As a result, increasing mg R (top-right frame) reduces \Qq\ and \Q' 6 \, while affecting 
IQ4I and \Q' A \ very little; recall that the latter two quartic charges receive the dominant 
contribution from the interference of t— and u— channel exchange diagrams. Increasing 
\fj,\ (bottom-left frame) has the same effect, as expected from our earlier observation that 
Q 6 and Q' 6 need sizable gaugino-higgsino mixing, while Q 4 and Q' 4 do not. Finally, the 
bottom-right frame shows that the dependence on the beam energy is relatively mild. 

Another conclusion from Fig. 2 is that \Q' 6 \ is usually the smallest of the four CP-odd 
quartic charges. The reason is that in this case t— and u— channel diagrams tend to 
cancel, whereas they add up in \Q$\. This indicates that measuring the production-level 
asymmetry will be quite challenging, as will be discussed in the next Subsection. 

6.2 Production asymmetries 

The simplest choice for probing the CP-odd contribution from Q' 6 to the production cross 
section in Eq. (17) is [8] 

/prod = sign[cos 9 sin(2$)] . (34) 
Instead a partly optimized asymmetry is suggested from the choice 

ip° r P od = cos6 sin 2 0sin(2$), (35) 

where we have set the angle 77 = for simplicity; nothing is gained by considering non- 
vanishing angles between the transverse e + and e~ polarization vectors. The factors of 
sin 2 and sin (2$) appear explicitly in the differential cross section in Eq. (17); inclusion 
of the factor cosG, which strictly speaking violates the construction principle described 
in Sec. 5, is necessary in this case, since this contribution to the cross section changes sign 
when cos G — > — cos G. 

Here it is appropriate to show that the asymmetries defined in Eqs. (30), (34) and 
(35) are indeed CP-odd. This can most easily be seen by using the so-called naive or T 
transformation, which inverts the signs of all three-momenta and spins, but (unlike a true 
T-transformation) does not exchange initial and final state. In the absence of absorptive 
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phases* a violation of T invariance is equivalent to CP violation, as long as CPT is 
conserved (which is certainly the case in the MSSM). Recall that we fixed the +z and +x 
directions via the e~ beam and spin directions, respectively, which are themselves T odd 
quantities.* In this coordinate frame a T transformation therefore amounts to flipping 
the signs of only the ?/— components of all three-momenta and spins. This is equivalent to 
flipping the sign of the azimuthal angle $ (as well as that of which is however irrelevant 
for the production-level asymmetry), leaving (and 0*) unchanged. Our production- 
level asymmetries are therefore T odd, which probe CP-violation if absorptive phases can 
be ignored. 

The effective asymmetries resulting from Eqs. (34) and (35) are shown by the (green) 
dotted and (black) solid curves, respectively, in three frames in Fig. 3. In these figures we 
have chosen the same default parameters as in Fig. 2, which ensures that x® ~* X\Z is the 
only possible two-body decay of As noted in Sec. 3, in this case we can measure the 
X® polarization only if the polarization of the Z boson is determined. In particular, one 
has to be able to distinguish between the two transverse polarization states in order to 
construct CP-odd asymmetries involving the Z polarization. Although this measurement 
is, in principle, possible for Z — > decays, the efficiency is quite low due to its small 
branching ratio (~ 7% after summing over e and [i final states), and a very poor analyzing 
power (from almost purely axial vector coupling for Z£ + £~). Although qq final states have 
larger analyzing power, the measurement of the charge is very difficult. It may be only 
possible to probe the production level asymmetry through this decay mode. 

Unfortunately the event cannot be reconstructed in this mode, as noted in Sec. 4. This 
means that we do not know the angles and $ appearing in the definitions of Eqs. (34) 
and (35); the best we can do is to approximate them by the corresponding angles of the Z 
boson. This leads to the (blue) dashed curves in the frames of Fig. 3 that show effective 
asymmetries, which are based on the "optimized" choice in Eq. (35). 

The top-left frame shows these asymmetries as functions of the CP-odd phase $i. 
We see that the "optimized" effective asymmetry exceeds the simple asymmetry based on 
Eq. (34) by typically ~ 20%, leading to a ~ 40% reduction of the luminosity required to 
establish the existence of a non-vanishing asymmetry at a given confidence level. Unfor- 
tunately replacing the true production angles (0 and $) by those of the Z boson reduces 
the effective asymmetry by a factor of 2.5—3.5. This suppression factor depends on the 
masses of the two lightest neutralinos, which in turn depend on $i. In this case even for 
the most favorable choice of parameters an integrated luminosity of several ab _1 would be 
needed to establish a non-vanishing optimized asymmetry at the la level, even assuming 
100% beam polarization! This is well beyond the currently expected performance of the 
international linear collider. 

^In the present context absorptive phases can only come from the finite width in the Z— propagator, 
which is entirely negligible for s 3> m|, or from loop corrections. 

■'■Note that for r\ = the initial state is T self-conjugate in this coordinate frame. 

§The effective asymmetry constructed from \% ~^ Xi^ decays is very similar to that from y^—*y^[Z 
decays; we therefore do not show numerical results for this decay mode. 
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Figure 3: The top-left and both bottom frames show the effective production-level asym- 
metries defined by Eq. (34) (green dotted curves, labeled "prod.") and (35) (solid black 
curves, labeled "opt. prod. "), together with the "optimized" production asymmetry where 
the true production angles are replaced by those reconstructed from the Z direction (blue 
long-dashed curves: without cuts; red short-dashed curve: with the cuts described in the 
text). The top-right frame shows the total cross section for e + e~ — > X1X2 without (black 
solid curve) and with (blue dashed curve) cuts. The default parameters are as in Fig. 2, 
but one parameter is varied in each frame. 

The lower-left frame of Fig. 3 shows that the situation might be better at higher 
beam energies. The effective production asymmetries peak at yfs ~ 900 GeV for the given 
choice of SUSY parameters. Moreover, the difference between the "theoretical" optimized 
asymmetry and the one constructed from the Z boson angles becomes much smaller at 
higher energy. The reason is that at ^/s ^> m^o the X2 becomes ultra-relativistic; its 
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decay products then fall in a narrow cone around the x° direction, so that the differences 
between the real production angles (G and $) and the corresponding angles derived from 
the flight direction of the Z boson become small. However, even in this case 1 ab -1 would 
only allow to establish an asymmetry with a significance of 3.5 standard deviations at 
best, ignoring experimental resolutions and efficiencies, and assuming 100% transverse 
beam polarization. The bottom-right frame shows that the situation is even worse if the 
mass of the SU(2) singlet selectron e R is close to that of the SU(2) doublet e^, which is 
taken as 500 GeV in this figure. 

The top-right figure is a reminder that X1X2 production can nevertheless provide useful 
information on the phase $1 [4], simply through a measurement of the total production 
cross section, which increases by almost a factor of three when $1 is varied from to 71; 
no beam polarization is needed for this measurement. As explained in Refs. [11, 4] this 
is due to the fact that the production occurs in a pure P— wave for $! = 0, but has a 
large S— wave component for $! = ir. This figure also shows that, for the chosen set of 
parameters, cutting against the ZZ background as described in Sec. 4, as well as applying 
the acceptance cut 

|cos6 x |<0.9 (36) 

for all visible final state particles X (in this case, the Z boson), only reduces the cross 
section by ~ 15%. The (red) short-dashed curve in the bottom-left frame shows that 
these cuts affect the effective asymmetries even less. 

6.3 Decay asymmetries 

We now turn to the "decay" asymmetries, which are sensitive to the x° polarization. We 
saw in Sec. 3 that these can be only probed through X2 ~ > Z^Z^ decays (ignoring three- 
body decays, which will be highly suppressed if any two-body decay is allowed). The 
discussion of Sec. 4 showed that in this case we can reconstruct the event with two- or 
four-fold ambiguity. 

Equation (21) shows that there are three CP-odd terms in the x\ polarization vector, 
which are sensitive to transverse beam polarization. In order to construct the correspond- 
ing "optimized" asymmetries, we first need an explicit expression for the scalar product 
appearing in Eq. (29). Working in the reference frame where the +x direction is defined 
by the transverse part of the e~ polarization vector, and using the same set of axes for 
the definition of the x° decay angles 0*, $* in the x° res t frame, we find using Eqs. (24) 
and (25): 

"ft-kl = V T [cos 6 sin6* cos($ - $*) - sin 6 sin 6*] 
+ V L [sin 6 sin 6* cos($ - $*) + cos 6 cos 6*] 
+ Pjvsine*sin($-<r), (37) 

where we have suppressed the superscript 2 on the components of the x% polarization 
vector. This, together with Eqs. (21) and (23), leads to the following choices for / in 
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Eq. (30):* 



f LN = [sin sin 0* cos($ - $*) + cos cos 0*] sin(2$) sin 2 , 

f TN = [cos 6 sin 0*cos($ - $*) -sin 6 sin 6*] sin(2$) sin(20) , 

f NT = [sin 0* sin($ - $*)] cos(2$) sin 9. (38) 

In each of the three expressions the factor in square brackets comes from Eq. (37), the 
second factor from Eq. (21), and the last factor from the expressions for E^, Stat an d 
Etvt, respectively, in Eq. (23). 

Similarly, the expression for A 2 ^ in Eq. (21) contains two CP-odd terms that can be 
probed with only longitudinal beam polarization, or even with unpolarized beams. Since 
the expressions for T^nu and Xtvl in Eqs. (22) and (23) are identical except for differ- 
ent quartic charges, we can combine these two terms into the "optimized" longitudinal 
effective asymmetry Al = A[fi] with 

f L = [sinB* sin($ - $*)] sin 9. (39) 

Note that the four functions fi defined in Eqs. (38) and (39) are all orthogonal to each 
other, i.e., the product of any two different functions will vanish when integrated over the 
entire phase space. 

Although the three asymmetries defined in Eqs. (38) are independent of each other 
(probing different £ab)> i n the context of the MSSM they all probe the same quartic 
charge Qq. If m^o and m^o are known, one can therefore construct a single asymmetry to 
probe Qq, called the total "optimized" transverse decay asymmetry At = A[/t] with 

f T = [sin O sin 9* cos($ - $*) + cos cos 6*] sin(2$) sin 2 ■ (l + ^ 2 - // 2 ) 
+ [cos sin 0* cos($ - $*) - sin sin 0*] sin(2$) sin(20) • /x 2 
+ [sin0*sin($ - $*)]cos(2$)sin0 ■ 2^ 2 , (40) 

where the /i, have been defined in Eq. (5). The first, second and third line in Eq. (40) 
correspond to the contributions from S^at, £ TA t and Etvt, respectively. 

Finally, we also consider an effective asymmetry based on the measurement of the 
momentum of the positive lepton l\ coming from the first stage of decay, defined by 
At = A[f+] with 

A + = sin(2$^). (41) 

The advantage of this asymmetry, which is somewhat similar to the decay asymmetry con- 
sidered in Ref . [8] , is that it does not need event reconstruction, as long as the "primary" 
and "secondary" leptons can be distinguished. 

As discussed in the previous Subsection, a CP-odd observable changes sign when 
$ — > — $ and $* — > — $*. Evidently the asymmetries defined in Eqs. (38) through (41) 

*Note that the denominator A^ 1 in Eq. (20) cancels against the factor A^ 1 from the production cross 
section (17) in the final result for the cross section differential in production and decay angles. 
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satisfy this condition. Due to the sign flip in Eq. (29) all asymmetries discussed in this 
Subsection have opposite signs for x° — ¥ £r&~ an d X2 ~^ ^r^ + decays; events of these two 
kinds should be treated separately. Since there are equal number of events from these two 
decay chains, we can simply focus on events with only positively charged primary leptons. 




Figure 4: Effective transverse decay asymmetries for the same default parameters as in 
Fig. 2, except that now rrie R = 155 GeV. The (black) dot-dashed, (magenta) long dashed 
and (blue) short dashed curves show the "optimized" asymmetries based on frN, Jnt an d 
/ln in Eq. (38), respectively, while the (red) solid curves show At of Eq. (40), and the 
(green) dotted curves show A± of Eq. (41)- In the right (left) frame acceptance and 
background-removing cuts have (not) been applied. 

The two figures in Fig. 4 show the effective "optimized" decay asymmetries based on 
Eqs. (38), (40) and (41). We use the same default parameters as in Figs. 2 and 3, except 
that the e# mass has been reduced to 155 GeV, so that X2 — * ^ e=F decays are allowed 
and dominant. Our choice of mg R implies that m^o — mg H ^> ni£ R — m^o. As discussed 
in Sec. 4 this implies that the harder lepton always comes from the first step of X2 decay, 
allowing to reconstruct the event with only a two-fold ambiguity. We average over both 
of these solutions when calculating the "optimized" asymmetries. We find that the wrong 
reconstruction typically leads to asymmetries with the same sign as the true solution, 
with (of course) smaller magnitude. The dilution of the asymmetries due to the event 
reconstruction ambiguity is therefore not very severe. The effective asymmetry based on 
$ln of Eq. (38) and, especially, the one based on f T of Eq. (40) are therefore substantially 
larger in magnitude than the simple effective asymmetry based on Eq. (41). Note also that 
the three effective asymmetries based on Eq. (38) move "in step" , as expected from our 
earlier observation that they all probe the same quartic charge Qq. Combining them into 
a single effective asymmetry, as in Eq. (40), therefore increases the size of the asymmetry 
significantly. 

The two frames in Fig. 4 differ in that the left figure does not include any cuts whereas 
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in the right figure we remove events that can be reconstructed as W or e~R pair background 
events. Also, we apply the acceptance cut in Eq. (36) to both final state leptons. For 
the case at hand these cuts only reduce the effective asymmetries by 10% to 20%. This 
high cut efficiency is also due to our choice of masses, which implies that the two leptons 
in the final state have very different energies. In contrast, both background processes 
have identical energy distributions for the two leptons in the final state. Signal events 
can be rarely reconstructed as background in this scenario. As a result we find that even 
after cuts one would only need an integrated luminosity of ~ 10 fb _1 to measure a non- 
vanishing asymmetry at the 3a level. This still assumes 100% beam polarization. Even 
for the more realistic choice PtPt — 0.5 one might achieve 3<r significance with ~ 40 
fb" 1 of data. This integrated luminosity should be achievable, assuming that transverse 
beams will be available. 

Finally, the four figures in Fig. 5 compare the simple asymmetry A± of Eq. (41), the 
total optimized transverse decay asymmetry At, and the optimized longitudinal decay 
asymmetry Al- We note that the longitudinal decay asymmetry is usually bigger than 
our total optimized transverse asymmetry. At least for probing the CP- violating phase in 
the context of the MSSM (where <3?i is the only relevant phase in the convention where 
M.2 is real), therefore, one does not really seem to gain anything by transverse beam 
polarization. The only exception is at large energy (bottom-right frame); this is due to 
the extra factor m^o/y/s appearing in the expressions for T, NU in Eq. (22), and T, NL in 

Eq. (23), which determine the size of Al. 

The upper right panel shows a quite complicated dependence of the effective asym- 
metries on me R . For intermediate cr masses both final-state leptons in signal events can 
have similar energies. As a result one often has four solutions for the event reconstruc- 
tion. In this case one cannot identify the "primary" lepton used in Eq. (41). We have 
dealt with this by simply discarding events with four solutions, since averaging over all 
four solutions would dilute the asymmetries a lot. Unfortunately this reduces the cross 
section significantly. At the same time e,R pair events become more similar to our X1X2 
events, since, as we just mentioned, the signal now has similar distributions for both final 
£ energies. Hence the cut against selectron pair production removes more signal events 
in the present case. As a result, the complete set of cuts reduces the total cross section 
by up to a factor of 5, the worst case being m e - fl ~ 195 GeV. Note that the different 
asymmetries are not equally sensitive to these cuts. The total "optimized" transverse 
decay asymmetry At is reduced by at worst a factor of 2, whereas the simple asymmetry 
A\ can go down by a factor of 4. The reason for this is that the cut efficiency depends on 
the same production and decay angles that appear in the definitions of our asymmetries. 

The lower left panel includes the longitudinal decay asymmetry Al for two different 
choices of longitudinal beam polarization. In both cases we take opposite polarization 
for the e + and e~ beams, since we are dealing with chiral couplings, see Eq.(ll). Usually 
taking a right-handed electron beam is most advantageous, since it maximizes the e.R 
exchange contribution; note that the e,R coupling to Binos, which is needed to probe the 
CP-odd phase $1, is two times larger than that of e~L- However, for very large this 
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Figure 5: Comparison of the simple transverse decay asymmetry (41) (green dotted 
curves), the total "optimized" transverse decay asymmetry (40) (red solid curves), and 
the "optimized" longitudinal decay asymmetry (39), the latter both for transverse (black 
dot-dashed) and for longitudinal (blue dashed) beam polarization. The default values of 
the parameters are as in Fig. 4, but one parameter is varied in each panel. 

choice is no longer optimal. In this case x° becomes more and more wino-like, i.e., it does 
not couple to e^. A right-handed e~ beam means that exchange does not contribute; 
the Z-exchange contribution also vanishes for large However, taking left-handed 
electrons one still gets a sizable contribution from ei exchange to the cross section, and 
also to the asymmetry. In the opposite regime of rather small \fj,\ the asymmetries depend 
very strongly on this parameter, since here X2 changes from a higgsino-like to a wino-like 
state. 

As in the previous figures (as well as in Ref. [8]) we took e^ beam polarizations ±1. 
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In the case of longitudinal beams one can then suppress the W or e# pair background 
(but not both), by appropriate choice of polarization. However, in practice the beam 
polarization will be significantly smaller than this; we therefore left the cuts against both 
backgrounds in place. We also note that longitudinal beam polarization can increase Al 
significantly, although the very small size of this effective asymmetry for our "default" 
parameters and transversely polarized beams (top left frame) is clearly accidental. 

Last but not least, we have checked numerically the effect of varying the left-handed 
selectron mass rrie L on the CP-odd asymmetries. The transverse decay asymmetries, with 
transversely polarized beams, are sensitive to the mass; in fact, they get a bit bigger with 
smaller mass values. Nevertheless, we have noted that the longitudinal asymmetry for 
unpolarized beams becomes much bigger when the left-handed selectron mass is reduced. 
For example, taking parameters as in the top-left frame in Fig. 5, except for a reduced 
mg L = 250 GeV, the maximal value of \ At\ after cuts increases to about 1.2 fb -1 / 2 , whereas 
the maximum of \ Ai\ reaches about 2.2 fb -1 / 2 . We emphasize that we do not actually need 
any beam polarization to probe this asymmetry, although it can be increased significantly 
by using longitudinal polarized beams; for reduced bl mass, taking left-handed e~ and 
right-handed e + beams is often optimal. Therefore, reducing the left-handed selectron 
mass does not affect the ordering of At and Al, i.e. the inequality Al > At (for optimized 
choice of longitudinal beam polarization.) 

7 Summary and Conclusions 

In this paper we studied the production of neutralino pairs at future linear e + e~ colliders, 
with subsequent two-body decays of the heavier neutralinos. We found that decays of 
the type Xi ~~ * X%^, Z) are not sensitive to the Xi polarization, unless one can measure 
the polarization of the Z— boson (or that of the final-state neutralino x*j)- These decays 
can therefore only be used to probe CP violation in neutralino production. Unfortunately 
the corresponding CP-odd term suffers from cancelations between t— and u— channel 
diagrams, and is nonzero only in the presence of higgsino-gaugino mixing. As a result, 
measuring this asymmetry, which can be done only with transversely polarized e ± beams, 
will be very difficult, if not impossible, with the currently foreseen linear collider perfor- 
mance. 

In contrast, xl decays into a slepton plus a lepton allows to probe the xl polariza- 
tion state, thereby opening up the possibility to construct several decay asymmetries. 
Moreover, this decay, followed by subsequent £ — > dxi decays, allows to reconstruct even 
the simplest neutralino pair events, X2X1 production with invisible (e.g., stable) x°, with 
two- or four-fold ambiguity. Under favorable circumstances experiments at a collider 
with (sufficiently strongly) transversely polarized beams should then be able to determine 
non-vanishing asymmetries with high statistical significance. However, even in this case 
a different asymmetry, which does not depend on transverse beam polarization (but can 
be maximized using longitudinal beam polarization), is generally larger in size than even 
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the best of the transverse decay asymmetries we studied. We saw in Fig. 5 that this is 
true both for gaugino- and higgsino-like It a l so remains true when we vary the ratio 
\Mi\/M2, in particular for \M\\ > Mi- However, if \M\\ ^> M2, or if both produced 
neutralinos are higgsino-like, all CP-odd asymmetries become small. Recall that in the 
MSSM all these asymmetries essentially result from a single (potentially large) phase, 
associated with the U(l) gaugino mass (in the convention where the SU(2) gaugino mass 
is real and positive). 

We therefore conclude that, at least in the context of neutralino production in the 
MSSM, transverse beam polarization is not particularly useful in probing explicit CP 
violation. Once the relevant masses have been determined, the most sensitive probe of 
the relevant CP-odd phases remains the total cross section [4], although it is a CP-even 
observable. If this measurement indicates that some phase differs from or ir, one needs 
to see explicit CP violation, in order to convince oneself that the variation of the cross 
section is indeed due to a phase, rather than due to some extension of the MSSM. However, 
as noted above, this can be most easily accomplished by using longitudinal, rather than 
transverse, beam polarization. 

The situation might be different in extensions of the MSSM, however. Whenever the 
quartic charges Qq and Q' 6 defined in Sec. 2.2 contain (combinations of) phases that are 
independent of those in Q4 and Q' 4 , the option of transverse beam polarization might be 
very useful for determining these phases. In the NMSSM, for example, the neutralino mass 
matrix contains additional CP-odd phases associated with the singlino sector, which can 
be large. A dedicated analysis along the lines presented in this paper would be required to 
decide whether transverse beam polarization could be helpful in disentangling this more 
complicated neutralino sector. 
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